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Abstract 

The emergent spatial patterns generated by growing bacterial colonies have been the focus of intense study in 
physics during the last twenty years. Both experimental and theoretical investigations have made possible a clear 
qualitative picture of the different structures that such colonies can exhibit, depending on the medium on which 
they are growing. However, there are relatively few quantitative descriptions of these patterns. In this paper, 
we use a mechanistically detailed simulation framework to measure the scaling exponents associated with the 
advancing fronts of bacterial colonies on hard agar substrata, aiming to discern the universality class to which the 
system belongs. We show that the universal behavior exhibited by the colonies can be much richer than previously 
reported, and we propose the possibility of up to four different sub-phases within the medium-to-high nutrient 
concentration regime. We hypothesize that the quenched disorder that characterizes one of these sub-phases is 
an emergent property of the growth and division of bacteria competing for limited space and nutrients. 

Introduction 

The use of statistical physics to unravel biological problems has yielded many novel insights into the living 
world. Numerous examples can be cited, including the application of statistical mechanics to ecological or 
metabolic networks [TJ [3] , or the use of phase transition theory and correlation functions to understand 
the behavior of animal collectives [3J 2]. One of the most remarkable instances of marriage between 
physics and biology concerns the emergent spatial patterns produced by growing bacterial colonies. 

Beginning with a strip or drop of bacterial culture inoculated on an agar plate, the growth and division 
of cells over time results in an emergent colony pattern that varies with nutrient and agar concentration 
in the underlying medium. This system presents many different biological and physical questions whose 
answers require the application of different conceptual approaches. For instance, we can analyze the 
influence of cell group structure on the evolution of cooperative or competitive secretion phenotypes, 
measuring the reproductive fitness of idealized strains [5j [6] . On the other hand, if we are interested in 
a morphological characterization of the colony's propagating front, different physical properties of the 
system must be analyzed. 

The Morphological Phase Diagram of Bacterial Colonies 

More than two decades ago, Matsuyama et al. [7] and Fujikawa and Matsushita [5] showed that bacterial 
colony patterns obtained in the laboratory can be fractal objects [9]. The properties of the emergent 
fractal depend mainly on two factors: i) nutrient concentration, C, strongly influences cell growth rate, 
and ii) agar concentration, C'a, controls cell mobility by altering the hardness of the substratum on which 
the colony is expanding. In the absence of special forms of collective cell motility [10] , bacterial colony 
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Fi gure 1; Sketch of the morphological phase diagram shown by bacterial spatial patterns (adapted with permission from 
|15| ; copyrighted by the Physical Society of Japan). Our work will be focused on the shaded zone, called here the rough 
regime. 



patterns may be classified with a two-dimensional phase diagram. Using Bacillus subtilis, Fujikawa, 
Matsushita and collaborators [HI [HUH] defined a phase diagram that comprises five regions with different 
qualitative behaviors resulting from the interaction between biological and physical factors (see Fig[T|): 

- At low agar concentration, cells can exhibit active flagellar motility. At low nutrient concentration 
(zone I), individuals tend to deplete local nutrient supplies rapidly and disperse in search of zones 
with better growth conditions. The resulting colonies consist of many long, thin branches (the 
so-called dense branched morphology, DBM) 13 . When nutrient concentration is increased, cells 
grow to higher densities with increased resource availability along the advancing front; the resulting 
pattern is compact and smoothly circular (zone II). 

- At medium agar concentration and relatively high nutrient availability, colonies consist of concentric 
rings that result from periodic alternation between static and motile stages with positive growth 
rates (zone III) [HQS]. 

- At high agar concentration, the substratum becomes hard and dry enough that cells cannot move 
by active means. Thus, low nutrient concentration results in patterns similar to the DBM, but 
with fewer, thicker branches due to limited cell movement jTTJ [T2J [T3] . These patterns (zone IV) 
resemble those created by diffusion- limited aggregation (DLA) processes [TB], with a very similar 
fractal dimension [7J[5]. For medium-to-high nutrient levels, the colony patterns are compact again, 
but with more irregular fronts than those obtained with low agar concentration (zone V) [121 E] ■ 

Many different models have managed to reproduce parts or all of the morphological phase diagram 
by adding new biological ingredients, such as cooperativity by means of chemotactic signaling [13] . or 
new physical ingredients, such as instabilities induced by fluctuations during cell diffusion [18) . Mat- 
sushita et al. 115] introduced a three-equation reaction-diffusion model capable of reproducing, at least 
phenomenologically, all regions of the phase diagram described above. 

Thanks to this considerable amount of both experimental and theoretical work, the qualitative picture 
of the bacterial colony pattern phase diagram is well understood. 

Quantitative Description 

Bacterial colony patterns are excellent candidates for analysis via propagating-front theory, due to a 
correspondence between the biology of bacterial growth and the physical aspects of nonliving fronts 
studied in statistical physics. Let us consider high agar concentrations (that is, increased surface hardness 
and minimized bacterial mobility) and an initially flat inoculum. As bacteria grow and divide, they 
passively shove each other, which causes the colony front to advance in space. Let h(x, t) be the distance, 
at time t, between the top of the cell group and the basal substratum at its linear coordinate x (see 
Fig©- Importantly, not all the cells can grow and divide. Only those that can access nutrients, namely 
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those forming an active layer at the advancing front, can do so (see inset in FigJ^J). The depth of this 
active layer depends on nutrient penetration into the group. The identification of growth and shoving 
as the main driving forces allows us to draw an analogy with non-living propagating fronts or driving 
interfaces from statistical physics, like those formed by fire on a burning sheet of paper or water absorbed 
into a porous medium [9]. 

The main observable measured in the study of these fronts is the variance in surface height, also 
termed the interface width or roughness (a measurement of the lateral correlations at the interface). It 
is defined as: 



W(L,t) = (\h(x,t) -h\ 2 \ , (1) 

where the overline represents spatial averages, < . > the average over runs of the experiment, and L is 
the inoculation strip length (system size, hereafter). A family of different 2-dimensional colony patterns 
may arise, depending on the environmental conditions. If the colony pattern is fractal or compact, the 
front can be treated a priori as a self affine object, and the so-called Family- Vicsek scaling may be 
assumed [20 : 

{t for t < t x 
(2) 
L a for t > t x 

That is, for short times the colony roughness follows a power law behavior defined by the growth exponent, 
j3; at time t = t x , it saturates to a size-dependent value W — W s t(L) whose finite-size scaling defines the 
roughness exponent, a. From Eq.@ we can deduce a power- law dependence of the saturation time on 
size: 

t x ~ i Q//3 , (3) 

where z = a /ft is the so-called dynamic exponent. Thus, two independent exponents determine the 
behavior of the interface. 

Vicsek et al. [T7J performed the first experimental measurement of one of these scaling exponents. 
They obtained a non-trivial value for the experimental roughness exponent, a exp , associated with the 
growth of Escherichia coli and B. subtilis. Some time later, similar values were independently obtained 
using B. subtilis in the rough regime 12 . 

In the experimental measurements reported in the literature described above, the a exp exponent was 
measured by choosing different parts of a growing colony as replicates at successive time points of 
their saturation regimes. At each of these time points, several boxes of lateral size I were used for the 
measurement of finite size effects over the roughness of the interface [T2l EZ1 ■ Let us briefly analyze this 
observable, which we designate w exp . If we define the "local roughness" of the front as [21] : 

w(l,t) = (<[&(*, t) - </^] 2 >*) 1/2 , (4) 

where < . >i represents a spatial average over boxes of size I, the behavior of this local observable 
resembles that of the global variable, but with different exponents: 

( t? fori<<4 
w(l,t)~{ . (5) 

[ for t > t' x 

Systems in which a ^ ai oc show what is called anomalous roughening, while regular self-affine interfaces 
fulfill the standard scaling for the roughness a = ai oc [35]. For / ~ a, where a is the grid spacing (or, 
in this case, the cell diameter), w(l ~ a, t) exhibits power-law behavior with positive exponent for cases 
of anomalous roughening, while it is constant or decreases with time for standard self-affine interfaces 
[2"Tj . For these standard self-affine interfaces, the finite-size scaling of this local roughness provides a 
good measurement of the global, universal exponent a. For interfaces showing anomalous roughening, 
however, the values of local exponents are not necessarily universal [23j . and we therefore cannot infer 
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Table 1: Value of the exponents for the different universality classes, experiments and numerical simula- 
tions mentioned in the text. The values with an asterisk have been calculated using Eq.([3]) or by using 
the standard self-affinity feature of w(l,t). The two values for a exp for the moving qKPZ correspond to 
measurements using large L (left) and a system size equivalent to those used in our simulation framework 
(right). 

the universal behavior of the system only by measuring them. In both cases, as I — > L, w(l, t) converges 
to W(L,t). 

The observable measured in experiments, w exp , is the result of performing this local finite-size scaling 
by averaging over multiple time points of the same run, using different parts of the same system as 
independent replicates. Given enough separation between time and spatial sampling points, w exp — 
w(l,t) and, thus, 

&exp — &ioc- In the experimental work mentioned above, measurements of this w exp 
performed in zone V yielded a scale-free behavior of the observable, with an exponent in the range 
&exp = 0.74 — 0.78 with a standard error of up to 0.07 [T^l [T7]. This result is supported by several 
theoretical models [2H[5S] (see table 1). 

Thus, with simple concepts and observables commonly used in statistical physics, one may character- 
ize the non-trivial morphology of bacterial colony fronts and quantitatively compare the colony patterns 
produced by different strains of bacteria. Similar values for the exponents measured in different bacterial 
strains indicate that, despite the varying microscopic details and interactions specific to each strain, they 
share the same basic biological and physical ingredients at a larger scale of observation. That is, they 
show the same universal behavior. 

Much less attention has been paid to the universal behavior exhibited by bacterial colony patterns. 
Attending to the expected symmetries in the colony patterns (i.e. the different terms of the equations in 
the theoretical models), the current consensus is that the branches of the low nutrient regimes are the 
result of a DLA process (zone IV) or DBM (zone I). The concentric-ring pattern (zone III) is treated 
as a special case, and the flat circular pattern (zone II) as a trivial one [25]. However, the behavior of 
the compact pattern observed in phase V (the rough regime, from now on) remains contentious. This is 
mainly because the exponent associated with the only measured observable, the roughness exponent, is 
insufficient to be conclusive in its classification. With no other observable measured, we cannot be sure 
if, for instance, the interface shows anomalous roughening, in which case a exp would not represent the 
global roughness of the front. In that case, the measurement of the local exponent would not provide 
any information about the universal behavior of the surface (see above). 

Here, we will attempt to fill this gap by studying bacterial colony fronts from the point of view of 
universality. By focusing on the rough regime (zone V) , we will show that the behavior of the propagating 
surfaces can be much richer than that described in the cited experimental work, prompting a modification 
of the phase diagram portrayed above. Also, we will provide the whole set of exponents defined in Eqs.([2])- 
(|3"|) supporting our classification, and determine which are the relevant biological and physical ingredients 
giving rise to such values for the exponents. 

A Proxy for the Experiments 

To measure experimentally the proposed set of exponents, we must be able to reproduce the environmen- 
tal conditions belonging to the regime under consideration, namely the rough regime, for each replicate. 
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Figure 2: Patterns resulting from the growth of the colony starting from an initial inoculum of L = 500 fim with periodic 
boundary conditions. Nutrient concentrations are C = 0.05 g/l (left) and C = 3 g/l (right). Inset: Illustration of the 
active layer (green circles) and how it distorts the bulk liquid nutrient field, creating a boundary layer; inactive cells are 
represented in black; reproduced from [5]. 



Here, cells are not actively motile, so growth and passive shoving are the sole mechanisms that induce 
colony front propagation. 

To characterize the behavior of the bacterial colony fronts, an ideal experiment should use independent 
replicates to measure the exponents presented in Eq. ([2])-Eq. (J5]) . An exhaustive sweep of the nutrient 
concentration values with sufficient replication to obtain good statistics is highly time-consuming, and 
though we are performing such experimental work (to be reported elsewhere |29|). here we will use a 
mechanistically detailed simulation framework. Our framework is derived from the latest generation of 
agent-based models developed by chemical engineers to predict the structure and metabolic activity of 
bacterial communities [30]. The model has already been used to study colony structure patterns in the 
high agar concentration regime (zones IV and V) , with excellent experimental support [301 123 122] ■ With 
this simulation framework, we can monitor the state of every cell and its interaction with the nutrient 
field and the rest of the group. Also, we have total control over environmental conditions and individual 
traits, including maximum growth rate and the ability to absorb nutrients |33j . The tractability offered 
by these individual-based simulations allows us to perform a finer exploration of the space of parameters 
that contribute to colony pattern formation. This will be useful not only for the goals of this paper, but 
also for guiding our ongoing and future work with laboratory experiments. 

Briefly, our framework comprises two overlapping layers, one of which is divided into a grid that tracks 
local nutrient concentration, and the other of which monitors bacterial cells that arc implemented off- 
lattice, i.e. as rigid circles in a continuous space. We consider the growth of a generic bacterial strain 
whose physiological traits are parametrized to those of E. coli. Cells grow according to local nutrient 
concentration, and once they reach a maximum radius of 0.75 /im, they divide into two daughter cells. 
A shoving algorithm is used to prevent overlapping between cells as they grow and divide. The model 
assumes a bulk fluid in which nutrient concentration is held constant, and implements a boundary layer 
in which nutrient concentration profiles are calculated by solving the following reaction-diffusion equation 
to equilibrium at each time step: 



where C/ oc is the local nutrient concentration, D is the nutrient diffusivity, Y is a metabolic yield coeffi- 
cient, and fj, is a Monod-type bacterial growth rate expression [33]. In calculating nutrient concentration 
profiles this way, we follow the common assumption that reaction-diffusion is much faster than bacterial 
growth and metabolism. More details on the simulation framework can be found in [30] . 

Fixing all other parameters, we will vary the bulk nutrient concentration, C, and system size, L, to 
measure the value of the exponents defined in Eqs.([2])-([S]). The system size will range from L = 300 fim 
to L — 600 /im, covering the maximum window sizes used in experimental measurements (up to 400 fim 
in [12j). In addition, we will monitor w(l ~ a, t), to discern whether the system shows anomalous 
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Figure 3: Left: Plot of the w eX p(l) function for two different system sizes; after an initial transient, they both converge 
to a power law with exponent a exp = 0.67(5). The colony pattern is compact and rough (inset). Right: Dynamic behavior 
of the W(i) function, showing a power law with an exponent = 0.61(5) for more than lOh of colony growth. Right, inset: 
Behavior of w(t); the zero slope of the curve indicates that a = a; oc = a exp . 



roughening (see above). Due to technical reasons, the duration of the experiment is limited to the time 
needed by the colony to reach h max = L, which depends mainly on nutrient concentration. 

We begin by checking that our model generates a qualitatively similar behavior to that found in the 
high-agar concentration experiments. Figure [2] shows that low (left) and medium (right) bulk nutrient 
concentration result in colony patterns similar to those of zones IV and V, respectively, with thick sparse 
DL A- like branches in the former case and a compact rough pattern in the latter HZ] • 

Now, we focus our attention on a quantitative characterization of the rough regime. For one selected 
nutrient concentration, C, belonging to this regime, we first measured the exponent a exp . One example 
of the w exp curve can be seen in the left panel of Fig|3] (red points), where an initial strip of L = 400 fim 
wide was grown with nutrient concentration C — 3.8 g/l (the resulting colony pattern is shown in the 
inset of the figure). Choosing this size, a saturating regime was found for the duration of the simulation, 
at which we measured a exp = 0.67(5). For the larger I, as can be expected, w exp saturates due to the 
approach of w(l, t) to W(L, t). This result is very robust, for it is shared by the whole range of compact, 
rough patterns tested (a exp = 0.62(5) — 0.75(5)). The values of the exponents measured fall into or close 
to the interval reported in experiments, within error bars (see table,). 

We also measured, for L — 400 /im, the local roughening exponent, obtaining ai oca i — 0.68(5). The 
temporal evolution of w(l <C L,t) (see right inset, Fig|3]) shows a constant behavior, indicating that 
ot = OLi oc — a exp for the colony. The growth exponent was measured for the largest system size (L = 
600 /im), giving a value of /3 = 0.61(5). Using Eq.Q, the derived dynamic exponent is z — 1.11(17). 
By collapsing the curves for different system sizes ranging from L = 300 fim to L = 600 fim [34], we 
independently confirmed these measurements. These exponents are in accordance with those found in 
[2"5] with a reaction-diffusion model similar to that introduced in [TO] (see above). We now turn our 
attention to determining the universality class to which the system belongs. 

A Universality Class for the Compact, Rough Patterns 

The Eden model is the paradigm of toy models representing growing cell colonies [35, 36, 37]. In one of its 
versions, a site along the colony front is randomly chosen each time step and one new cell is placed at any 
of its free surrounding locations [36 . The basic ingredients of this model are i) a source of stochasticity 
working each time step, i.e. an annealed noise, due to the random choice of the site at which to deposit 
a new cell, and ii) lateral growth, generated by the possibility of choosing a neighboring site for the 
new cell. These two main ingredients are described by the mesoscopic Langevin equation known as the 
Kardar-Parisi-Zhang equation [26] : 

d t h{x, t) = X(\7h) 2 + D\7 2 h + F + ar)(x, t), (7) 
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where D, A, F and a are constants, and 77, an uncorrelated (white) noise. The first term in the r.h.s. 
represents lateral growth and the second term represents an effective diffusion generated by it. The 
exponents characterizing this universality class have been calculated not only numerically, but also 
analytically; in 1 + 1 (x and h) dimensions, they are: a = ai oc — 1/2, ft — 1/3, z — 3/2 [9]. 

Due to the similarities between the Eden model patterns and bacterial colonies, KPZ has been suggested 
as the universality class that describes colony growth. Most of the theoretical models for bacterial growth 
require a non-linear diffusivity function for the active cells to generate compact patterns similar to those 
of phase V. The simplest of these functions, in which diffusion is proportional to the local density of 
active cells, develops eventually the lateral growth and effective diffusion associated with the growth of 
the front in the KPZ equation, Eq.0 [55] . 

However, the large a exponent found in experiments does not match that of the KPZ class (see table 
1). It represents stronger lateral correlations, and has been argued to be the result of nonlocal interaction 
between the moving cells, or a consequence of the rod shape of the bacteria used in experiments (E. coli 
and B. subtilis). In the rough regime, these particular strains form folding chains lying very close together, 
causing lateral correlations to develop more quickly [12] , Such a large a exponent value is also found in 
systems with different stochastic terms, like non-Gaussian thermal noises, or quenched noises [9l 139] . 

If we replace the thermal noise of Eq.© by a quenched noise, T}(x, h), the interface changes its universal 
behavior. For F < F c , the advance of the surface is hindered and eventually prevented by pinning 
sites, represented by the quenched noise. For F > F c , the interface is able to overcome the pinning 
places and grows indefinitely. At F = F c , the interface undergoes a second order phase transition 
between the pinned (non-moving) phase and the depinned (moving) phase, characterized by the exponents 
a = a ioc = 0.63(3), ft = 0.67(5), z = 1.01(1) @D] (or a = 3/4 ft = 3/5, z = 5/4, analytically deduced 
in [27]), which define the so-called quenched KPZ class (qKPZ). This scaling is retained in the moving 
phase (a = ai oc = 0.75(5), ft = 0.62(3), z = 1.18(5), own measurements). Note the large value for a, 
close to those observed in the experiments and numerical simulations of the theoretical models (see table 
Although different sources of quenchedness such as heterogeneities of the agar substratum [5J |2~4] 
have been proposed to influence expanding bacterial colonies, the exact origin for the large a exponent 
is still controversial. 

Interestingly, our mechanistic model produces a roughness exponent similar to those reported in real 
experiments, but none of the proposed explanations for this exponent are present: the cells are circular, 
shoving events move them randomly very short distances, and there is no possibility for substratum 
heterogeneity in our framework. Moreover, our measurements agree with the set of theoretical exponents 
for the qKPZ class. Cell-cell shoving events generate the short-range lateral growth and diffusion terms 
and the thermal noise but, what is the source of quenchedness? 

In our case, quenchedness is an emergent feature of the colony. It is the result of a non-trivial interaction 
between cells and the nutrient field gradient. The growth of each individual is entirely dependent on 
the amount of nutrient that it can access (see above), which is conditional on the interaction between 
the nutrient field and the focal cell and its neighbors. Individuals residing at the peak of a region in the 
colony with larger height than its immediate surroundings have better access to nutrients and grow more 
quickly. On the other hand, growth rate is reduced among individuals at the bottom of regions with 
smaller height than the surroundings, where nutrient concentration is reduced due to consumption by 
neighboring cells; the propagation of the colony front is therefore reduced in these "valleys" (see Fig0]). 

Competition between cells for the available resource can thus generate points at the interface with 
substantially reduced propagation velocity |41j . Moreover, this inequality of propagation is subject to 
positive feedback: the slower the propagation of the front at a site in comparison with the neighbors, 
the deeper the valley becomes and the more severely nutrient access decreases. This feedback effect can 
cause individuals in the valley to stop growing entirely, creating effective pinning sites along the front. 
Such pinning sites can only be overcome by the interface by means of an increased external pulling force 
(nutrient concentration, C, in this case) or by lateral growth, which are the same mechanisms exhibited 
by qKPZ interfaces. 

For the values of bulk nutrient concentration used for our measurements (see above) , the pinning sites 
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Figure 4: Generation of effective pinning sites by means of the feedback mechanism described in the text. The nutrient 
concentration at the bottom of the valleys is so low that cells stop growing. Cells at the top, however, can access the bulk 
concentration, growing at a much larger rate (see 1331 V Left: State of the colony at time t = 41h. Right: state of the 
colony 8 hours later. 



are always eventually overcome, behavior that corresponds to the moving phase of a propagating front 
advancing through a random environment (quenched disorder) . In this way, we can see that the bacterial 
profiles of our agent-based model have the same basic ingredients and quantitative behavior of a qKPZ 
interface in its supercritical phase. For system sizes similar to those of our framework, the theoretical 
finite-size exponents of that class agree with our measurements (see Table 1). 



A Richer Rough Phase 

In addition to the quenched noise, we must consider the stochasticity inherent in the shoving events 
among growing and dividing cells for a full description of cell colony patterns. We will aim to explain 
how the competition between these two noises could act as the nutrient concentration is varied in the 
rough regime, by using the general equation: 

d t h(x,t) = \{Vh) 2 + DV 2 h + F + a im (x,t) + a 2 r] 2 (x 7 h), (8) 
with F and the ratio o~\jo~i changing with C: 

- For a medium nutrient concentration, the resource is quickly depicted in valleys, while the rest of 
the front can grow at different velocities. Pinning sites therefore emerge easily, while the number 
of shoving events is low because the growth rate and, thus, the number of birth events, depend on 
nutrient concentration |33j . In Eq.([5]), this is translated into a small o~\jo~i- Therefore, this source 
of annealed noise remains asymptotically irrelevant for the macroscopic behavior of the front as 
compared with the self-generated quenched noise. It is the reported moving-qKPZ regime. 

- As nutrient concentration is increased, the growth rate of the cells also increases and so does 
the mean velocity of the interface. Birth events and the subsequent shoving are more frequent. 
This hinders the formation of deep valleys, eventually overcome by the advancing colony front. 
This corresponds to both a larger F and a larger amount of annealed noise (larger a\), leading 
to a behavior that would be governed by the KPZ equation due to the irrelevance of the present 
quenchedness. 

- If the bulk nutrient concentration is further increased, most or all cells in the colony layer have 
access to inexhaustible resources and grow at the maximum rate at all times (F> 1), leading to a 
completely flat colony surface. A summary sketch can be found in Figj5] 



These observations motivate further exploration of the parameter space, varying the value of the 
nutrient concentration in search for the different regimes described above. While we could find a range 
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Figure 5: Sketch of the rough part (zone V) of the phase diagram, showing the different sub-phases discussed in the text. 



of concentrations with exponents in agreement with those of qKPZ universality class (see Table 1), we 
were not able to measure, for any C, a set of exponents fully compatible with those of the KPZ class. This 
might be due to the time limitation of the framework: the higher the value of bulk nutrient concentration, 
the shorter the time required by the growing colony to reach h max (see above), making it more difficult 
to look for crossovers in high C regimes. We may, thus, attribute the missed KPZ regime to a lack of 
observation of the asymptotic behavior of the observables. We are currently developing a new version of 
the simulation framework in order to overcome this technical limitation [29] . 

At very high nutrient concentrations, a completely flat interface is easily recovered. The transition from 
a flat behavior of W(L,t) to a non-trivial value for its associated exponents indicates that the bacterial 
colony undergoes a roughening phase transition. However, such flat examples will be hardy seen under 
real experimental conditions (see below). 

At very low bulk nutrient concentration the interface cannot overcome pinning sites, while the rest 
of the active layer continues growing. This behavior is shown in FigLU it does not correspond to the 
supercritical phase of a qKPZ system, nor its subcritical phase, because the interface never becomes 
completely inactive. For such low nutrient concentrations, the equivalence with surfaces propagating in 
random environments is not valid anymore. In this regime, the large differences in velocity between the 
growing parts of the front and the pinned sites make the /3 exponent take unusually large values around 
j3 = 1 (with patterns similar to those reported for the "groovy" phase in 42 ). At even smaller bulk 
nutrient concentrations, we enter a phase in which wide pinning regions are developed, but the colony 
continues to grow at a slow pace. Here, branches arise and the colony can grow not only upwards and 
laterally, but also downwards due to the presence of more ample free space. This description corresponds 
to patterns that belong to zone IV in FigfT] 

Conclusions 

Bacterial colonies are among the best examples of biological systems whose complexity can be informa- 
tively reduced by using concepts from physics. Here, we have shown that the essential biological and 
physical ingredients shaping colony propagating fronts may be distinguished by quantifying the scaling 
behavior of a colony. 

By using a detailed individual-based model in lieu of real experiments, we have shown that cell shape 
and substratum heterogeneities — in contrast to the existing literature — are not necessary to explain 
the form and behavior of colony fronts in the rough regime. We have also hypothesized that this regime 
contains a much richer morphological variety than expected. Based on the biological and physical 
ingredients present in our simulation framework, we have given biological justifications for the emergence 
of quenched disorder as the result of competition for space and resources among cells. We have further 
described how higher nutrient concentration allows colonies to overcome pinning sites due to higher cell 
growth rates, giving rise to different universal behaviors. Overcoming the technical limitations of the 
simulation framework will grant better access to the asymptotic behavior of the front in each of these 
sub-phases. 

Classifying the universal behavior of a generic bacterial colony is more than a mere academic exercise. 
It allows us to describe the behavior of the front by means of one single equation, such as Eq.©, instead 
of using explicit agent-based or reaction-diffusion theoretical models. We have observed that irregular 
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cell shape, long-range cell motility, and extracellular compounds are not necessary for our idealized cell 
groups to exhibit the same universal behavior as real bacterial colonies. This result suggests that such 
factors may not be responsible for the behavior of the colony patterns observed in experiments, and that 
the same universal behavior may be expected for the morphology of any bacterial colony under similar 
conditions. 

However, for the moment, we do not have enough information to prove this claim. Only experimental 
work with different strains will allow us to check if the richer universal behavior described here can 
be observed even in the presence of other ingredients that can be potentially important for the real 
asymptotic behavior of colony fronts |29j . For instance, due to the average mutation rate of many 
commonly studied bacteria, mutations that affect growth rate often arise on the time scale reported for 
colony growth experiments that we have discussed. A rare minority of these mutants exhibits higher 
growth rates than the rest of the population, resulting in bursting sectors that emerge from growing 
colonies. These sectors obviously influence the shape of colony fronts in a way that is not considered by 
any theoretical models in the literature so far. These mutations will prevent also the observation of the 
roughening transition in real experiments: the larger the value of C , the larger the number of replication 
events and the higher the likelihood of obtaining mutants with a larger growth rate. This could generate, 
by different means, large lateral correlations leading to a large a exp such as that measured in the cited 
experimental studies. 

All of the open questions discussed here highlight the necessity of a complete quantitative charac- 
terization of the propagating of real bacterial colonies in the rough regime. With it, the biology of the 
interaction between individuals and the physics of the environment, studied from a multidisciplinary 
point of view, will allow us to discern the relevant mechanisms responsible for the asymptotic behavior 
of real bacterial colony profiles. 
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